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Abstract. We use simple models (the Ising model in one and two dimensions, 
and the spherical model in arbitrary dimension) to put to the test some recent 
ideas on the slow dynamics of nonequilibrium systems. In this review the focus 
is on the temporal evolution of two-time quantities and on the violation of the 
fluctuation-dissipation theorem, with special emphasis given to nonequilibrium 
critical dynamics. 



Prologue 

The aim of this review is to summarise recent works devoted to the dynamics of 
ferromagnetic spin systems after a quench from infinite temperature to their critical 
temperature. 

The initial impetus for such an investigation was the desire to put to the test, 
on simple models, some recent ideas on the slow dynamics of nonequilibrium systems 
(aging of two-time quantities and violation of the fluctuation-dissipation theorem). 
By simple models we mean models with no quenched disorder, with, for some of them 
at least, the virtue of being solvable. Here we address the case of ferromagnetic spin 
systems, such as the Ising model in one and two dimensions, and the spherical model 
in arbitrary dimension. Urn models are also simple enough to serve the same purpose. 
They are the subject of another review in this volume llj). 

During the course of this investigation we realised the interest of posing the same 
questions for nonequilibrium critical dynamics [^J ||. 

1. The fluctuation-dissipation theorem and its violation 

Consider a generic spin system evolving at constant temperature from a disordered 
initial configuration. 

Let s and t, with s < t, be two successive instants of time, and r = t — s, their 
difference. Denoting by a(t) the spin at time t, we consider the correlation 

C(t,s) = (a(s)*(t)), 
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and the local response to a time-dependent external magnetic field H (t) 

At equilibrium, that is when the waiting time s is large compared to the 
equilibration time r oq , these functions are stationary. They only depend on the time 
difference r: 

C(s,t) = C eq (r), 
R(t,s) =i? oq (r), 

and are related by the fluctuation-dissipation theorem (for a simple presentation see 
e.g. |): 

This situation is typical of the high-temperature regime (e.g. T > T c for a 
ferromagnet), where r cq is small. 

In experiments or simulations, instead of measuring R(t,s), one considers the 
integrated response, i.e., either the thermoremanent magnetisation of the system at 
time t, -Mtrm(^ s), obtained after applying a small magnetic field h, constant between 
t = and s; or the zero-field-cooled magnetisation Mzpc(ij s), where now h is constant 
between s and t. Defining the reduced integrated response p(t, s) by 

P (t,s) = -^M(t,s), 

we thus have 

PTKu{t,s)=T duR(t,u), (1.1) 
Jo 

Pzfc(M) —T J duR(t,u). 

At equilibrium, using the fluctuation-dissipation theorem, we have 

,C(t) 

Ptrm(M)= / dC = C(r), 



PzFc(t,s) = [ dC = l-C(j), 
Jc(t) 

thus a plot of p against C is given by a straight line of slope +1 (/>trm) or —1 (/9zfc)j 
as soon as s is large enough. 

At low temperature (below T c for a ferromagnet), r cq is either very large or 
infinite. In the scaling regime where 1 <C s ~ t <C T cq , aging takes place, i.e., C and i? 
are no longer stationary, and the fluctuation-dissipation theorem does not hold. The 
question is therefore to determine the relationship between C and R, if any. This can 
be done by defining the fluctuation-dissipation ratio X(t, s) by ||, |[ [?J 

^, S ) = — 

Assume that, in the scaling regime, all the time dependence of R can be 
parameterised by C. Or, in other words, that C acts as a clock for R. That is, 
for 1< s ~ t, 

A(M)~A(C(M))- (1.2) 
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As a consequence, we have 



,C(M) 

Ptrm(M)« / dCX(C), 
Jo 

PZFc(M) « / dCX(C). 



'C(t,s) 

Hence, in a plot of p against C, the slope at a given point is given by ±X(C). 

This behaviour has been observed in a number of instances. In particular, a census 
of the different cases of spin systems hitherto studied shows the existence of three main 
types of behaviour at low temperature (for a summary, see ^] , and references therein) . 
For domain-growth models, X(C) is discontinuous in C, taking a first value equal to 1, 
and a second one equal to zero fl Q (see discussion in section 2). For spin-glass 
models with p-spin interactions, X(C) is still discontinuous but the second value is 
non-zero. Finally, for continuous spin-glass models, X{C) is a non-trivial curve ||. 

In the present review we show that, at T = T C1 non-trivial statements can be 
formulated on the same issue. Hereafter we specialise to ferromagnetic spin systems. 
We take as representatives the Ising model in one and two dimensions, and the 
spherical model in arbitrary dimension. The Hamiltonian describing these models 
reads 

where the first sum runs over pairs of neighbouring sites. 

For the Ising model, <Ji — ±1, and the (non-conserved Glauber) dynamics is 
governed by the heat-bath rule: 

V(<Ti(t + dt) = ±1) = - (1 ± tanh/9fci(t)) , 

where the local field reads hi = Y]^ <jj + Hi, the sum running over the neighbours of 
site i. 

For the spherical model, <jj is a real number with the constraint o~f = N, 
where TV is the number of spins [l^, [l4| . The dynamics is governed by the Langevin 
equation jn| 

do - .; BE . , . , . 

In the right side, X(t) is a Lagrange multiplier ensuring the constraint, and Tji(t) is a 
Gaussian white noise with correlation (r]i(i)r]j (t 1 )) — SijS(t — t'). 

In both cases, at time t — 0, the system is in a disordered initial configuration 
(e.g. corresponding to equilibrium at infinite temperature). 

2. Aging below T c : low-temperature coarsening 

We first describe in more detail the behaviour of correlation and response at low 
temperature, for a generic ferromagnetic model such as the spherical model or the 2D 
Ising model, evolving at constant temperature after a quench from T = oo down to 
T <T C . We defer the discussion of the ID Ising model to section 4. 

In such a situation, domains of opposite sign grow, with a characteristic size 
L(t) ~ i 1 / 2 , where z — 2 is the growth exponent 17|. 
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In a first regime (1 ~ r < s), dynamics is stationary. Correlations decay from 
C(s, s) = 1, to the plateau value 

<7ea = lim lim C(s + t,s) = M| , 

r — >oo s — >oc H 

where M e<1 is the equilibrium magnetisation. Though the system becomes stationary, 
it is still coarsening, and therefore does not reach thermal equilibrium. However the 
fluctuation-dissipation theorem holds, and X = 1. 

In the scaling regime where s and t are simultaneously large (1 < s ~ f), with 
arbitrary ratio x — t/s, aging takes place, and correlations behave as |l7|] 

C(M)«M c 2 q /c(£). (2.1) 

For small temporal separations (r <C s, or x — > 1), we have fc(%) ~ ¥ lj implying 
C(t,s) — > M^ q . In other words, equation (2T) describes the departure from the 
plateau value M 2 q . For well-separated times (1 < s < t, or i > 1) fc{%) decays 
algebraically as 

fc(x)^A c x-^ z , 

where A is the autocorrelation exponent ]T^ | . As a consequence, we have 

dC(t,s) _ Ml 
ds s 

with fc(x) ~ Ac x~ x l z , at large x. 

In the same regime it is reasonable to make the scaling assumption (see discussion 
below) 

R(t, s )^s- 1 - a f R (-) 1 (2.2) 



■fc>[- 



with an unknown exponent a > 0, and with again the decay at large x 

f R (x)^A R x~ x / z . (2.3) 
We have therefore 

l ' >~ M^ 1 f c ,(t/s)~ M|/ Ac,' 

The fluctuation-dissipation ratio thus vanishes in the scaling regime, irrespective of 
the ratio t/s. 

For instance, for the spherical model, the equilibrium magnetisation reads 

M c 2 q - 1 - J 
and the correlation C(t, s) is given by (|]l]) with 

fc{x) = (^ttf) ' 

hence the autocorrelation exponent A = D/2. The response is given, in the scaling 
regime, by 

R(t,s)^^)- D ^r-\ (t-s)- D /*, (2.4) 
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which is in agreement with the form ( |2.2| ), with scaling function 

f R {x) = {^)- D ^x D /\x-l)- D '\ 

and the exponent a = D/2 — 1. 

For the 2D Ising model, the exponent A « 1.25 [|l8| is only known numerically. 
This is also the case of the scaling functions fc and fn [^0|. The latter work is 
compatible with a = 1/2, as predicted in Refs. JTTj, plj ], where it is argued that the 
integrated response scales as p(t,s) ~ L(s)~ 1 g (L(t)/L(s)) for soft spin models with 
non-conserved dynamics. 

In summary, for a ferromagnetic spin system |l(| [H| [n| [ll| , 

• for short times (r <C s), such that C(t, s) > M^ q , the fluctuation-dissipation theorem 

holds, and X = 1; 

• for long times (r ~ s), such that C(i, s) < M| q , the fluctuation-dissipation theorem 

does not hold, and X(t, s) — > independently of the ratio i/s. 

Note that we have 

=S(C- M e 2 q ) , 

in agreement with the static interpretation of X(C) in terms of the distribution of 
overlaps P(q) ||. 

3. Aging at T c : critical coarsening 

The system is now quenched from T = oo down to T c . 

In such circumstances, spatial correlations develop in the system, just as in the 
critical state, but only over a length scale which grows like t l / Zc , where z c is the 
dynamic critical exponent. On scales smaller than t x l Zc the system looks critical, 
while on larger scales the system is still disordered. For instance, the equal-time 
correlation function C r (t) = (oo(i)ay (t)) scales as 

where (3 and v are the usual static exponents. (A summary of the values of exponents 
is given in the Table.) The scaling function g(y) goes to a constant as y — ► 0, while it 
falls off very rapidly when y — * oo. 

The same temporal regimes, as defined in the previous section, are to be 
considered. However, their physical interpretation is slightly different, since the order 
parameter M| q vanishes and symmetry between the phases is restored. 

In the first regime (r <C s), the system again becomes stationary, so that the 
fluctuation-dissipation holds. 

In the scaling regime (r *~ s), temporal correlations behave as^ 

C(M)^" c /c(j), a c = 2[3/uz c = {D-2 + r,)/z c . (3.1) 

It is instructive to relate this behaviour to that observed for T < T c , namely, 
C(t,s) « Mg q /c (V s )- The passage from one formula to the other one is done by 

% For simplicity we use the same notation fc, fa, etc. for the scaling functions appearing in this 
section, though they are different from those appearing in the previous section. We use the same 
convention for the amplitudes Ac, An, etc. 
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noticing that in the critical region one has M eq ~ \T — T c \@ ~ £,cq^ u - Replacing £ eq 
by s 1 ^" implies the replacement of M c 2 q by s - 2 P/ uz c ^ s -(D-2+ri)/z c _ 

At large time separations (x 3> 1) we have (see |2j| for a derivation in the case 
of the so-called model A p4fl) 

where A c is the critical autocorrelation exponent ]2j| , related to the initial-slip critical 
exponent C p3| by A c = D — z c Q c . 
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Figure 1. Log-log plot of the critical autocorrelation function C(t, s) of the 
two-dimensional Ising model, against time ratio x = t/s, for several values of 
the waiting time s. Data are multiplied by s 2 /" z c . in order to demonstrate 



collapse into the scaling function fc(%) °f e( i- (3-iL. Straight line: exponent 
-Xc/z c « -0.73 of the fall-off at large x. (After ref - ] 0.) 



As a consequence of (3.1), we have 



dC(t, s) 
ds 



fc 



with the decay fc( x ) ~ Ac x x "l z " at large x. 

In the scaling regime, the response function behaves as 



R(t, s) 



,,, [ \ 



and, for large temporal separations, 
f R (x)KA R x- x °/*'. 



(3.2) 
(3.3) 



Nonequilibrium critical dynamics of ferromagnetic spin systems 



(Se e p3[ for a derivation of (B^) an d (|3.3| ) in the case of model A.) Note the similarity 
of (3.2) and (3.3) with (2.2) and (2.3), respectively. The scaling form (3.2) of the 
response implies 



pTRM(t,s) « s a " fp[ - 



with, as x > 1, f P (x) « A p a; -W z <\ 

We finally obtain the fluctuation-dissipation ratio 



X(t,s) 



T f R (t/s) =x (l 



and, at large temporal separations, 



A, 



lim lim X(t, s) = lim X(x) — T c 

s — >oo t — >oo x — >oo Aqi 



rp Ap 

C A C 



The last equality is equivalent to saying that, for 1 -C s -C t, 

PTRM(t, s) f» XooC{t, s). 

The limit fluctuation-dissipation ratio X^ can thus be measured as the slope near 
the origin of the C — ptrm plot. The scaling function X{x), and in particular the 
amplitude ratio X^,, are universal, in the sense that they neither depend on initial 
conditions nor on the details of the dynamics |j| . 

In the scaling regime, neither Ptrm(^ s) nor X(t, s) are functions of C(t, s). 
Instead, X(t,s) and s° c p(t, s) are functions of x = t/s, which is in contrast with 
the situations where equation (L2) holds, and further described in section 1. 

We now illustrate the results presented above. For the spherical model (see the 
Table for the value of exponents), the two-time correlation function reads 



where 



T, 



C(t,s) w s" 



Thus 



fc(x) = < 



K = 



T, 



c (D-2)(x + l) 
2(4ir)- D / 2 



(D/2 - 1] }c{x\ 

^^{x-l) 1 - ' 2 



D 



( K (x~l) 1 - D ^~(x + l) 1 - D ^ 2 ) 



2 < D < 4, 



D > 4. 



3D/2-2 2<£><4, 
D D > 4. 



Similarly, the response function behaves as 

R(t,s) : 



-D/2 



where the scaling function fii(x) reads 

(4^)- Z5 / 2 x 1 - Z5 / 4 (x - 1)- D ' 2 2 < £ < 4, 



Finally 



X, 



(4ir)- D / 2 (x-l)- D / 2 

1-2/D 2 < D < 4, 
1/2 D>4. 



D > 4. 
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In x 

Figure 2. Log-log plot of the critical integrated response function Ptrm (i-> s ) of 
the two-dimensional Ising model, against time ratio x = t/s, for several values 
of the waiting time s. Data are multiplied by s 2 P/ VZa , in order to demonstrate 
collapse into the scaling function f p (x). Straight line: exponent — X c /z c R3 —0.73 
of the fall-off at large x. (After ref. M.) 



For the 2D Ising model A c w 1.59 [g5| is only known numerically. Figures 1 
and 2 show numerical determinations of the scaling functions fc and f p Q . In two 
dimensions, we have X^ ss 0.26, and a preliminary study leads to Xoo w 0.40 in three 
dimensions [Q. 

The above discussion can be summarised as follows: 

• For short times (r <C s), such that C'(t, s) S - 2 P/ UZ " j the fluctuation-dissipation 

theorem holds, and X = 1. 

• For long times (r ~ s), such that C(t,s) ~ s ~2/3/vz c ^ fl uc tuation-dissipation 

theorem does not hold. The fluctuation-dissipation ratio X(t, s) is given by the 
scaling function X(t/s), such that X(x) — > X M as x — > oo. 

This is the critical counterpart of the behaviour of X(t,s) = X(C) for T < T c , 
summarised at the end of section 2. 

A last comment is in order. At thermal equilibrium, for a ferromagnetic system 
at criticality, the relationship between magnetic field and magnetisation, h ~ M^ q , is 
nonlinear. Therefore linear-response theory, used above to extract the response of the 
system, only holds for a magnetic field small compared to the scale ho ~ s -P&/vzc ^ 

s -(D+2- V )/2z c _ 
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4. One-dimensional Ising model at T = 

The one-dimensional Ising model is special in the sense that its critical temperature 
T c is zero. Hence the low-temperature phase does not exist. 

Another peculiarity of the model stems from the fact that the magnetisation 
exponent (3 is equal to zero. As a consequence, at criticality (i.e., at T — 0), there is 
no temporal prefactor in the expression of C(t, s) (or equivalently, no spatial prefactor 
in that of C r (t)). Indeed, let us recall that, at criticality, for a generic ferromagnetic 
model, we had 



C(t,s)^s- 2 ^f c 



t 

For the ID Ising model at zero temperature we have 

Irl 



C r (t) « erfc 



2 / 2s \ 5 

C(t, s) « - arctan ( - — -J . (4.1) 

The latter formulas are compatible with the former ones, taking into account that 
(3 = for the ID Ising model. Otherwise stated, the absence of an anomalous 
dimension implies that C(t, s) is not small in the critical region, in contrast to the 
generic cas es co nsidered in the previous section. 



From (4T) we obtain 
fc> (x) 



tt(x + 1) V x — 1 

The critical temperature T c being equal to zero, we define the dimensionless 
response function 

*«■•>-«■ 

In the scaling region (1 <C s ~ t), this function is found to behave as 

_! . . t 



R(t,s)^s- L f Ii [- 

where 



n^/2(x - 1) 

This again is compatible with the generic case, with (3 = 0. 

The reduced magnetisation ptrm(*, s) and the fluctuation-dissipation ratio 
X(t, s) can be computed explicitly. Both quantities only depend on t/s, or equivalently 
on C, in the scaling regime. One finds, in this regime ||, p6| , 

,~s V2 ( 1 ttC 

/Otrm(C ) = — arctan I —= tan — 

while A is more simply written in terms of the ratio x — t/s as 
Yd A - - x + 1 
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We note once again that the fact that (3 — implies no dependence in s in these 
quantities. Finally, the last equation implies for the limiting ratio 

Xoc = lim lim X(t, s) = -. 

S — ►QO t — >00 Z 



5. Discussion 



At criticality, for the generic cases of the spherical model and of the 2D Ising model, 
X{t,s) is not a function of C(t,s). It is instead a function of the ratio x = t/s, or 
equivalently of s 2t3 ^ Zc C(t, s) = fc{x). In this last representation, the value of X at 
the origin is equal to X^. Then the fluctuation-dissipation ratio increases and reaches 
the limit value 1 when the abscissa fc(x) goes to infinity, that is, for x — ► 1, where 
the fluctuation-dissipation theorem holds. 

Is the amplitude ratio I M related to equilibrium quantities? This remains an 
interesting open question. More generally, do the above results on the fluctuation- 
dissipation ratio admit a static interpretation, e.g. in terms of the distribution of 
overlaps P(q) [p2| ? Strictly speaking, the existence of a non-trivial X^ should imply 
the presence of an unexpected discrete component in P(q). We mention a recent work 
on related matters ^7j, where the finite-size behaviour of P(q) for the 2D X-Y model 
is related to the finite-time behaviour of p(t, s). 



A recent analysis 20 , based on conformal invariance, predicts the following form 
of the response function 

R(t,s) = r (t-s)- A (^\ , (5.1) 

without predicting the values of exponents appearing in the right side. This prediction 
should hold for a large class of systems. We note in particular that, for the spherical 



model, equations (2.2) and (3.2), together with the explicit forms of the scaling 
functions fn(x), given in sections 2 and 3, confirm this prediction, which is also verified 
by numerical computations on the Ising model in two and three dimensions |20[ . The 
analytical results for the ID Ising model given in section 4 do not, however, satisfy 



the prediction (5.1) 



Finally it is worth adding a few words on the comparison between the results 
reviewed here and those reviewed in ref. |l]] for urn models. For the zeta urn model, 
the situation at criticality is in all aspects similar to that of a generic ferromagnetic 



model, as described in section 3. However the prediction (5.1) is not fulfilled by this 
model. In the low-temperature phase, the results obtained for the zeta urn model do 
not fall in the framework reviewed in section 2, valid for a coarsening system. Finally, 
the results obtained for the backgammon model at T = arc rather different from the 
generic behaviour of a ferromagnetic model. A natural explanation of the discrepancy 
between urn models and ferromagnetic models is that in the former case the system 
is rather subject to condensation than to coarsening. 
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Table 1. Static and dynamical exponents of the ferromagnetic spherical model, 
and of the Ising model in one and two dimensions. First group: usual static 
critical exponents n, j3, and v (equilibrium). Second group: zero-temperature 
dynamical exponents z and A (coarsening below T c ). Third group: dynamic 
critical exponents z c , A c , and © c (nonequilibrium critical dynamics). 



exponent 


spherical (2 < D < 4) 


spherical (D > 4) 


2D Ising 


ID Ising 


T) 








1/4 


1 




1/2 


1/2 


1/8 





V 


l/(D-2) 


1/2 


1 


1 


z 


2 


2 


2 




A 


D/2 


D/2 


w 1.25 




Zc 


2 


2 


w 2.17 


2 


A c 


3D/2- 2 


D 


« 1.59 


1 


e c 


1 - D/A 





» 0.19 
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